The conductance of molecular wires and DFT based transport calculations 
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The experimental value for the zero bias conductance of organic molecules coupled by thiol- 
groups to gold electrodes tends to be much smaller than the theoretical result based on density 
functional theory (DFT) calculations, often by orders of magnitude. To address this puzzle we have 
analyzed the regime within which the approximations made in these calculations are valid. Our 
results suggest that a standard step in DFT based transport calculations, namely approximating the 
exchange-correlation potential in quasistatic nonequilibrium by its standard equilibrium expression, 
is not justified at weak coupling. We propose, that the breakdown of this approximation is the most 
important source for overestimating the width of the experimentally observed conductance peak and 
therefore also of the zero bias conductance. We present a numerical study on the conductance of 
the organic molecule that has recently been studied experimentally by Reichert et. al} that fully 
agrees with this conclusion. 

PACS numbers: 



I. INTRODUCTION 

Recently, several conductance measurements of single 
organic molecules have been reported^"^. Fig. 1 illus- 
trates a typical setup with one of the molecules used in 
experiments^. In part the tremendous attention that this 
field "molecular electronics" has received in recent years 
is due to its possible technological impact. Indeed, a 
molecular storage device has already been built^ and also 




FIG. 1. Schematic representation for a conductance mea- 
surement of the molecule (9,10-Bis((2'-para-mercaptophenyl) 
-ethinyl)-anthracene): between gold contacts [1]. 

molecular switches have been realized^. It is clear that for 
technological applications being able to systematically 
model and understand the charge transport properties 
of the combined system of molecule and contacts is of 
crucial importance. 

Since the conductance is very sensitive to details of 
spectral and orbital properties of molecules and their 
wavefunctions, ab initio methods like density functional 
theory (DFT) are an indispensable tool not only for 
structure calculations but also for transport theory. 
Nowadays, DFT calculations including several thousand 
electrons are possible, which allows the treatment of 
large molecules and also to include parts of the contacts 
for more realistic calculations. Despite these enormous 
capacities, a quantitative description of transport for 
weakly coupled molecules with a conductance well below 



the conductance quantum, .g<Cl, has not been achieved. 
In fact, experimental and theoretical values for the zero 
bias conductance of organic molecules, e. g. benzene, 
often differ by 1-3 orders of magnitude®'^. 

In the present work we consider the organic 
molecule depicted in Fig. 1, which has been sub- 
ject to experimental^ and DFT based theoretical 
investigations^", before. In this case the experimental 
value for the zero bias conductance is smaller than the 
theoretical one by a factor ^ 20. Three different reasons 
that possibly could lead to this discrepancy should be 
mentioned and discussed^ ^. 

First, Fig. 1 shows an idealized situation and ac- 
tual microscopic conditions realized in the experiment 
are not known. However, we have chosen the particular 
case. Fig. 1, because the experimental results are well 
reproducable^". If important variations in the atomic 
structure of the contact would exist, one should expect 
strong fluctuations in the conductivity from sample to 
sample. Indeed fluctuations are present, see Fig. 4 be- 
low, but their magnitude is relatively small so that strong 
structural variations are unlikely to occur. Small varia- 
tions in structure, on the other hand, tend to have little 
impact on the theoretical conductance, only. Our calcu- 
lations show, that the transmission is not very sensitive 
to changes in the bond angles or bond lengths as long 
as the change is within reasonable limits. Therefore, de- 
viations from the assumed atomic arrangement depicted 
in Fig. 1 from its experimental realization do not seem 
to offer a plausible explanation of the large discrepancies 
observed. 

Second, in theoretical calculations approximations 
have to be done when the molecule is coupled to the leads. 
Artefacts can be excluded only, when parts of the leads 
are included into the calculation ("extended molecule"), 
so that the following hierarchy of inequalities is met: 



<5gM < 7oM *C 7m 



(1) 
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{5cM' level spacing of extended molecule, 7m (Tcm) level 
broadening of bare (extended) molecule when coupled to 
the leads). Indeed, these conditions have not been met in 
previous calculations where i5eM ~ 7m- In order to im- 
prove upon this result, we have performed conductance 
calculations for the extended molecule. Fig. 1, where 
<5eM/7M ~ 0.1 and Ssm/Sm ~ 0.01 so that the hierarchy 
of inequalities is satisfied and artefacts from modeling of 
the leads can be excluded. We find that qualitative agree- 
ment of our data with the experiment is improved signif- 
icantly as compared to the earlier calculation^^. How- 
ever, the quantitative disagreement is not decreased but 
rather increased by another factor of 10! This shows that 
insufficient modeling can be ruled out as a source for the 
discrepancy. 

After we have discarded two obvious possibilities for 
explaining the observed discrepancy, we have to resort 
to more fundamental considerations. In this paper, we 
advocate a third possibility, namely that the standard 
implementation of DFT based transport calculations re- 
lies crucially on assiimptions that are not justified in the 
limit of "weak coupling" where the spatial structure of 
molecular wavefunctions is strongly inhomogenous and 
transport is dominated by single resonances. -"^^ 

We give a outline of our paper and a brief account 
of key results. In section II we discuss the derivation 
of the basic equation of DFT-transport, Eq. 11, and 
show that it has a wide range of applicability. This is 
true provided that the DFT-calculations are performed 
with an exchange-correlation potential appropriate for 
the nonequilibrium situation. However, since this poten- 
tial is not known in general, in practice an "equilibrium 
exchange correlation approximation" (EXCA) is made in 
which the standard equilibrium potential is used. Apri- 
ori, this approximation is uncontrolled. In section II B we 
introduce an alternative - but equivalent - formulation 
of transport, the Kubo formalism. It enables us to give 
analytic arguments for the applicability of the EXCA. 
Our results suggest, that the EXCA can be trusted only 
in the case of non-resonant transport, when the molec- 
ular level broadening is strong and individual molecular 
levels strongly overlap. A typical example for this case is 
a chain of gold atoms. By contrast, if the transmission 
exhibits resonances, transport is determined by individ- 
ual molecular orbitals and their broadening. In this case, 
the corrections to the EXCA become significant and in 
particular the level broadening and therefore the zero 
bias conductance may be severely overestimated. The 
molecule in Fig. 1 is a representative for these systems. 

In section III, we present numerical transport calcula- 
tions for a gold chain and the molecule, depicted in Fig. 
1, that we have already alluded to above. In view of 
our theoretical analysis we propose that it is the break- 
down of the EXCA in the weakly coupled limit that is 
the cause of the large discrepancy between theoretical 
and experimental molecular conductance. A discussion 
of our findings will be given in section IV 



II. TRANSPORT FORMALISMS 

A. Landauer-Buttiker-Formalism for Interacting 
Electrons 

Meir and Wingreen have derived a general expression 
for the current flowing through a region of space, where 
the charge carriers can interact - like a quantum dot or a 
molecule^^. Specifically, the Hamiltonian has the struc- 
ture 

H = Hoi{dl};{dm})+ J2 

a— /,r 

+ H.C.). (2) 

Tna=lr 

The first term describes the bare molecule which in gen- 
eral may include the electron-electron interaction. Its 
detailed structure will be of no importance in what fol- 
lows. The molecule makes contact to two leads denoted 
left and right. The bare leads are assumed to be nonin- 
teracting and described by the second term in Eq. (2). 
The third term represents the contact. 

The Meir- Wingreen formula connects the retarded (ad- 
vanced) Green's functions G (G^) and the lesser function 
of the full many-body problem (including the leads) 
with the dc-current 

1= jdE tT{hT^-fnTn){G - Gt) + tr(rL-rR)G< 

(3) 

where /l,r = ,f{E — /Ur^l) denotes the Fermi distribu- 
tion functions for the leads at chemical potentials /iR.L 
and Fl.r is the imaginary part of the self energy that 
describes the coupling of the molecule to the external 
leads. (For the most part the dependence on energy E 
will be suppressed in our notation.) For noninteracting 
particles, in terms of these self energies we have G~^ = 
EI-Ho-'Si^-'Er, sothatrL,R = «(Sl,r-SI„r) and 

G< = zG ( /l Tl + /r Tr ) Gt. (4) 

In the appendix, we demonstrate that Eq. (4) is just 
the statement, that the density matrix can be con- 
structed from left and right going scattering states V;,r 

{x = (x,f)), 

G<{x,x') = Y, fi M^) ri{x')+fr M^) rAx'). (5) 

l,r 

From (3) together with (4) we find for the transmission 
in the limit of linear response and zero temperature 

T{E) = tr TLGFRGt. (6) 

Eq. (6) is an incarnation of the familiar Landauer- 
Buttiker formula^* and valid for noninteracting electrons. 
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In this section, wc argue that expression (6) continues 
to hold also for a much larger class of interacting prob- 
lems. This is because according to the Runge-Gross the- 
orem of time dependent DFT (TDDFT) the time evolu- 
tion of any many-particle Hamiltonian can be calculated 
by solving a single particle problem in an appropriate ef- 
fective potential^^. Since this point is of importance to 
us we elaborate on it in the following. 

The effective time dependent single particle problem 
we should solve is of the form 

-l^^Onix) = i-^V' + V,{x))er^{x) (7) 

where 

Vs{t) = Vi + Vuln] + Vxc[n] + V.^t) (8) 

with Vi denoting the ion-core potential. Vh the Hartree- 
interaction, Vxc the exchange-correlation potential and 
14x the external probing field. (Notation suppresses the 
spatial index x.) In addition to the explicit time de- 
pendence of the potential Vs imposed by T4x an implicit 
dependence exists because Vg is a functional (in general 
nonlinear and nonlocal in time and space) of the den- 
sity n{x,t). Let us assume that Eq. (7) describes the 
molecule together with the leads and an external per- 
turbation V^x that is switched on at t = to and time 
independent thereafter: 

V,^ix) = v;x(x)0(t - to) 

At times prior to the system is assumed to be in 
thermal equilibrium (T = 0), so that the reservoirs are 
characterized by occupation numbers /„. Then we have 

'n{x) = Gv{x,x), with 

g<ix,x') = ^uen{x)e:ix') (9) 

n 

where 0„(x, f) evolves according to Eq. (7). More- 
over, we specialize to the case, where 14x(x) generates a 
monotonous electrical potential drop from V^x = V in the 
asymptotic region of the left lead to 14x = in the right 
lead. In response to the potential drop a current is being 
generated. After an initial period exhibiting transient be- 
havior, there will be a parametrically wide time interval 
in which the current and the density are quasistationary. 
This is precisely the situation for which also Eq. (3) has 
been derived^^. It is only at even much longer times, 
that the electrochemical potential becomes homogenous 
again and the current stops to flow. Formally speaking, 
we perform the order of limits in which the size of the 
reservoirs is send to infinity first and to — cxd, there- 
after. We mention that, since we are interested in the 
long-time-behavior only, details of how the external po- 
tential is switched on, step-like or adiabatic, are unim- 
portant. The corresponding memory is erased inside the 
reservoirs. 



Since Eq. (9) describes the exact evolution of the time 
dependent density we can also calculate the (longitudi- 
nal) current density j = dn/dt and hence find an exact 
expression for the conductance. Indeed, the same rea- 
soning that is used for noninteracting electrons to relate 
transport to a scattering problem can be employed for the 
present effective single particle problem as well. Conse- 
quently, a zero temperature description of the quasista- 
tionary region in terms of scattering states should exist 
with Gy taking a form similar to Eqs. (4) and (5) 

g<{x,x') = J2 ^'(^) ^tW) + fr erix) 9;{X'). (10) 

l,r 

It involves the Kohn-Sham orbitals 9r,i representing the 
scattering states of the quasistationary nonequilibrium 
situation and their zero temperature occupation numbers 
fi^r imposed by the left and right reservoirs. Since also 
the derivation of Eq. (4) can be repeated for an effec- 
tive single particle problem, a relation analogous to (3) 
also holds for Q'^ and the corresponding retarded Green's 
function Q. Together with the previous statement (10) 
this implies that the transmission is given by 

T = tr rLavrR^I^ (11) 

with Oy^{E) = Sl — Sr,, and the resolvent ma- 

trix 

where the sum is over KS-energies e„ and orbitals 0„(x) 
calculated for the uncoupled molecule. 

Eq. (11) constitutes the main result of this section. 
Similar arguments can also be found in a recent commu- 
nication by Stefanucci and Almbladh^^. 

Several aspects of our finding should no pass by with- 
out further comment, a) Eq. (11) constitutes the gener- 
alization of the Landauer formula to interacting electron 
systems. In the special case, where Fl and Fr, differ 
by a constant factor only, a derivation has already been 
given by Meir and Wingreen before^'^. We emphasize, 
however, that this condition is extremely restrictive. It 
implies that every atom of the molecule is coupled in pre- 
cisely the same way to the left lead as it is coupled to the 
right lead. Given the fact that physical couplings decay 
with increasing spatial distance, the "condition of pro- 
portional couplings" is violated for every realistic system 
with a finite extent. 

b) As before the self energies 5]l,r account for the cou- 
pling of the molecule to the leads. It is very well known 
that they can incorporate sophisticated many body ef- 
fects, like i.e. the Abrikosov-Suhl resonance if Kondo- 
physics is present However, one can restructure the 
problem and define an "extended molecule" that also 
comprises parts of the contacts, see Fig. 1. The new 
contact surface iSeM can be arranged sufficiently far away 
from the physical contact, so that the new self energy 
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depends on the type of lead only, but is totally inde- 
pendent of which molecule is used. In fact, if a separa- 
tion of energy scales exists such that the level spacing 
of the extended molecule, ^eivi, is much smaller than the 
broadening, 7m, of the bare molecular levels upon cou- 
pling to the leads the microscopic information carried by 
Sl,r becomes irrelevant. With any choice of 7cM in ac- 
cord with Erj. (1) the simple replacement I]l,r(x,x') 
= *7cM (^xx' with X, x' situated on iScm is justified. This 
freedom merely reflects the experimentalists choice to at- 
tach leads at convenience from any sort of shape or mate- 
rial as long as the voltage drop is near the real molecule. 
In this limit, all the detailed information about transport 
properties of the molecule including correlation effects is 
carried by the resolvent matrix (11) and inherited from 
the exchange correlation potential. 

c) Eq. 11 has been derived under the condition of van- 
ishing temperature and frequency and in the regime of 
linear response. Under these restrictions, scattering is 
energy conserving and therefore an effective single parti- 
cle scattering formulation of transport can exist. Upon 
releasing the constraints the incoming particle can ex- 
change energy with the molecule. Qualitatively new phe- 
nomena can occur and in general the scattering problem 
will become much more complicated and may not be un- 
derstood in terms of a simple single particle picture. For 
example memory effects appear, because here the incom- 
ing electron sees the molecule in the state it has been left 
in after interaction with the previous electrons. 

d) Eq. (12) suggests that near equilibrium (i.e. for the 
linear response) a variational principle may exist that al- 
lows for the calculation of the voltage at a given current 
or vice versa. The idea is to introduce a density operator 
that maximizes the entropy under the constraint that the 
current be finite. Work along this line has been done by 
Kosov^^, recently. 



B. Kubo Formalism of Linear Response 

The practical usefulness of Eq. (11) is limited since 
not much is known about the exchange correlation po- 
tential Vxc(^) that defines the Kohn-Sham problem 
in quasistatic nonequilibrium. Throughout all prcvioiis 
works using the DFT approach to transport it has been 
universally assumed that Vxc(^) may be approximated 
by its equilibrium form Vxc{0) used in standard DFT- 
calculations (EXCA). Let us discuss now, under which 
conditions this approximation may be expected to hold. 

Instead of solving a scattering problem one can also 
find the conductance from the alternative, but equiva- 
lent, Kubo formalism. The advantage of this starting 
point for our purposes is that analytical statements about 
the excitation frequencies are available which contain in- 
formation about the corrections to the EXCA. 

We observe, that the current is related to the dynam- 
ical polarization 



I = dP/dt. 



(13) 



As a consequence the linear current response to a ho- 
mogenous external electric field is governed by the dy- 
namical polarizability tensor of the entire system includ- 
ing molecule and leads 



I = —IU)CX{U)) Ee 



(14) 



which is closely related to the full density susceptibility 
aij{w) = J dx dx! Xi x'j x(x, x',a;) (15) 
with the corresponding Lehmann representation 



^0 ^-E^o+iO 



(16) 



The sum is over all many body states with energies 
starting from the ground state /i=0 and the abbreviation 
Ef_to=Ef^—Eo has been used. In the case of strongly reso- 
nant transport the sum is dominated by the contribution 
of a few poles and the off-resonant current results from 
the finite lifetime of the corresponding excitations which 
is encoded in the dipole matrix elements. Clearly, this is 
the regime in which molecules attached with thio-groups 
find themselves. 

The full density response contains a piece that de- 
scribes dynamical screening. In terms of the bare KS- 
response xks(x, x',a;) the relation 



X 



'\'^) = Xks(^) -vq- /xc(w) 



(17) 



holds true, where the exchange-correlation kernel 
fxc{x,x') = dVxG{x)/dn{x') has been introduced and 
i^Q = |x — x'|~^. In fact, Eq. (17) can be thought of as 
a definition of Vxc- The crucial point is that the addi- 
tional terms in (17) shift the excitation energies fl^ away 
from their bare KS-values = Ci — ej (we introduce a 
multi-index /U = (i,i)) 

% ^ a;f + (Ml VQlt,) + (/x|/xc(a;f)|M) (18) 

where = denotes a single particle transition from 

the occupied KS-orbit j to the vacant orbit i. When /xc 
is decomposed into its exchange contribution and the cor- 
relation part /c further analytical progress is possible^^ 
and one finds 



0„ 



,KS 



+ vq\h) - {ii\ VQ\jj) + (mI/cIm)- (19) 



,,KS| 



{v^^: exchange part of the bare Hartree-Fock (HE) op- 
erator, i;^^: exchange potential in the KS-equation). Eq. 
(19) contains the plausible statement that the actual ex- 
citation energies move away from their KS-values closer 
to their Hartree-Fock estimates when dynamical screen- 
ing is taken into account^^. Note, that a corresponding 
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change must be expected to happen also in the dipole ma- 
trix elements of Eq. (16) which can be important since 
they control the lifetime of the excitations. 

Equipped with this information we can proceed now 
and discuss the validity of the EXCA underlying previous 
DFT-transport calculations. Two limiting cases should 
be distinguished: 

i) individual resonances strongly overlap, 7m ^ Sm- 
This is the limit of strong coupling, g ^ 1, where the 
wavefunctions \i) are extended and the deviations from 
KS- and HF-excitation energies in Eq. (19) are small. In 
fact, for structureless, plane- wave-like states they vanish. 
In this case the standard approximation is justified. An 
example for a system in this class is a linear chain made 
from Au-atoms. 

ii) individual resonances are fully developed, 7m ^ i^m- 
Here, the wavefunctions show pronounced localized fea- 
tures and in general the corrections to KS-excitation en- 
ergies and matrix elements will not be small. While qual- 
itatively correct results may still be found in many cases, 
a systematic quantitative analysis of the resonant struc- 
tures in T{E) based on the EXCA is not possible, in gen- 
eral. A typical representative of weakly coupled systems 
is given in Fig. 1. The contact of the extended 7r-system 
to the leads is interrupted by an S-atom and therefore 
the molecular states are partly localized. 



III. NUMERICAL CALCULATIONS 



elements, tL.R(XA'^, xn) which describe a hopping pro- 
cess of an electron in an orbital state N of an atom at 
position X of the extended molecule to a state n of the 
atom at a location x in the leads: Sl — ^LflL^L 
similarly for Sr. The hopping matrix elements we ap- 
proximate by their bulk values that we obtain from an 
independent DFT calculation for a large gold cluster (146 
atoms). Likewise, we replace the surface Green's function 
(xn|gL,R|x'n') of the leads by a bulk one taken from the 
same calculation^^. 

B. Results 

Fig. 2 shows our results for the transmission of a linear 
chain of four equidistantly placed Au-atoms {d = 2.67^). 
As a check we have performed calculations with leads 
that have been modeled by 54 (see Fig. 1) and 84 Au- 
atoms . Moreover the number of "surface" Au-atoms 
from the extended molecule that have been coupled to 
external leads has been varied: the self energy has been 
calculated with 29 and 41 atoms taken into account. It 
can be seen that the transmission is (essentially) indepen- 
dent of these parameters, as it should. Our results are 
in very good agreement with experiments and previous 
calculations^'^. 

Before we present our results for the transmission func- 
tion the coupling of the organic molecule to the contacts 
should be discussed. 



The discussion presented in the previous section sug- 
gests that the experimental conductance of Au-chains is 
well described by the standard DFT-approach whereas 
the conductance of organic molecules coupled by thio- 
groups is only poorly captured. In the following section 
we present transport calculations making use of Eq. (12) 
together with the EXCA in order to corroborate this re- 
sult. 



A. Method 

Our method is similar to approaches described in^^~^^, 
however it has the advantage that using the program 
package TURBOMOLE we can include a considerably 
larger number of contact gold atoms. For the molecules 
of interest to us, this number (110) is sufficient so that 
the Fermi energy of the extended molecule is very close 
to the bulk value E-p — 5.1eV, even without attaching 
additional leads. 

Since the jest of the method has been outlined 
previously^ ^ ^"^ wc can limit ourselves to a brief descrip- 
tion. The transmission is given by Eq. (11) but with 
a resolvent matrix Qq that has been obtained employ- 
ing the EXCA: KS-orbitals and KS-cncrgics arc taken 
from a standard DFT calculation^^'^^. The self energies, 
Sl,r can be expressed in terms of the hopping matrix 
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FIG. 2. Transmission of a four atom gold chain. Upper 
panel: Contact was made of 54 gold atoms, see contacts in 
Fig.l. Out of these 41 atoms (□, 5 x 5 and 4x4 bottoms 
layers) have been coupled to external leads. For comparison, 
results with 29 coupling atoms are shown as well (o). Lower 
panel: Same as upper panel with 84 gold atoms defining the 
contacts. {E-p w —5.1) Traces are (nearly) independent of 
the modeling parameters chosen, and they agree well with 
experiments and previous calculations^'^ 
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FIG. 3. Transmission of the molecule in Fig. 1 with sulfur 
atom coupling to one (solid) , two (dashed) and three (dotted) 
gold atoms. {Ef ~ — 5.1eV) Modifications in the atomic con- 
act structure lead to a slight shift of the conductance peak 
that offer an explanation for the inhomogenous peak broad- 
ening observed in experiments, see Fig. 4. 

In experiments a standard way to facilitate this coupling 
is to introduce a sulfur atom that forms a very strong 
bond with gold as well as with carbon atoms. Since 
the precise microscopic conditions of experiments are not 
known, several possibilities for the sulfur-gold bonds have 
to be considered. It is known, that sulfur tends to bind 
to three gold atoms on a plane Au(lll) surface, i. e. 
the hollow site is the most stable one^*. Binding to just 
one gold atom - on top position - corresponds to a lo- 
cal minimum of the free energy and binding to two gold 
atoms is unstable. We find that the situation on rough 
surfaces exhibiting edges, is different. Near an edge, the 
sulfur finds its most stable position by binding to two 
gold atoms (Fig. 1) in agreement with an earlier study 
on much smaller systems^^'^°. The on top position re- 
mains a minimum albeit at higher energy (-|-0.7eV) and 
the hollow position is unstable. Because in the break 
junction experiments of interest to us the sulfur atom is 
likely to be exposed to an irregular surface we consider 
all three cases. 

In Figs. 3 and 4 we show the transmission and the I/V- 
characteristics of the molecule of Fig. 1. The traces cor- 
respond to the S-atom binding to one, two or three Au- 
atoms. These micro-deformations induce a slight shift of 
the transmission peak near — 5.2eV by about 0.2eV. 

Qualitatively, the observed features arc in accord with 
experimental findings of [1] , inset of Fig. 4: we find a con- 
ductance gap and the differential conductance exhibits 
a maximum at about 0.4V that stems from a resonant 
molecular level - the remnant of the HOMO of the bare 
molecule - about 0.2eV below the Fermi energy. The peak 
conductance and the zero bias conductance can differ by 
an order of magnitude and at larger voltages the conduc- 
tance rapidly decreases before there is another increase, 
again. On the quantitative level serious discrepancies 
between our calculations and experiments persist: our 
value for the zero bias conductance 0.2e^//i exceeds the 
experimental one, ~ 10~^e^//i, by more than two orders 
of magnitude. Moreover, the width of the experimental 
peak is strongly temperature dependent and can decrease 
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FIG. 4. Upper curves: dl/dV curve from data in Fig. 3. 
based on /(V) = JdET{E) {f{E + V/2) - f{E - V/2)) with 
f{E) denoting the Fermi-function (T = 300K). Lower curves: 
experimental data from [1]. Different traces represent re- 
sults from consecutive voltage sweeps. The theoretical and 
experimetal traces exhibt the same qualitative features, but 
strong quantitative deviations in the conductance magnitude 
exist. 

by an order of magnitude with temperature decreasing 
from 300K to SOK^o. This strongly suggests that the 
experimental peak width at room temperature is de- 
termined by inhomogeneous broadening rather than a 
strong coupling to the leads. 

The inhomogenous broadening can be understood on 
the basis of our results in Fig. 3. We have seen above 
that the position of the HOMO is roughly accounted for 
by our DFT procedure and it is mainly its damping that 
is overestimated. Therefore, the shift of the peak position 
upon changing the coupling of the molecule to the leads 
is indicative of a microscopic smearing mechanism: at 
sufficiently large temperatures the thermal average over 
the different types of S-Au-bonds leads to an effective 
broadening on the energy scale « 0.3eV which indeed is 
in accord with experimental observations. 

One might suspect that at least part of the reason for 
the large discrepancy between theoretical and experimen- 
tal findings for the zero bias conductance is that in real 
experiments the molecule is exposed to various boundary 
conditions, e. g. stress, that prohibit the contact geome- 
try to relax completely, which is what we have assumed 
in our calculation. Therefore, an important question is 
whether the conductance can be strongly affected by a 
slight change in microscopic degrees of freedom, like the 
S-Au-bond length or the bond angles defined in Fig. 5. 




FIG. 5. Schematic representation of geometric degrees of 
freedom defining the sulfur-gold-coupling, (center atom: sul- 
fur, left: gold surface, right: benzene-molecule) 
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FIG. 6. Transmission of benzene- 1,4-dithiol for various pa- 
rameters a, d and /3 as defined in Fig. 4. Plot shows, that the 
transmission is robust against small structural modifications 
induced e.g. by strain. Upper panel: molecule, coupling to 
a Au(lll) hollow site, fully geometry optimized (solid) and 
same molecule after rotation (dot-dashed, a — 7r/24). Mid- 
dle: change in bond length (optimal value: d = 2.55A). ap- 
plied change O.OSA (dashed) and 0.1 A (long dashed). Bottom: 
change in angle (3 by 7r/12. {Ef « — 5.1eV) 

Since our interest is in a qualitative question, we per- 
form our study using the simpler molecule benzene- 1,4- 
dithiol. Our findings for the transmission as a function of 
energy are in excellent agreement with recent results by 
Xue and Ratner^^ and Stokbro et al.^^. In Fig. 6, upper 
panel, we display the transmission of this molecule for the 
geometry optimized case and also after a subsequent rota- 
tion of the molecule about angle a — 7r/24 defined in Fig. 
4. The middle panel shows the impact of changing the 
sulfur-gold bond length d by 0.05A and O.lA, the lower 
channel exhibits the change upon changing the bond an- 
gle /3 by 7r/12.. These manipulations have only a small 
effect on the transmission confirming earlier findings^^'^^. 
The example shows that the DFT-conductancc docs not 
change by orders of magnitude when changing details of 
the sulfur-gold bonds within reasonable limits. 

Now, we arc facing the following situation: On the 
one hand, the present calculations show that small vari- 
ations in the nanostructure do not strongly influence the 
transmission. On the other hand, experiments are re- 
producible and do not show very strong fluctuations in 
the zero bias conductance, which implies that large vari- 
ations in the atomic structure of molecule and contact 
do not exist. Furthermore, the chemical bond between 
S and Au is known to be very strong and stable this 
is why sulfur has been chosen as the coupling element in 
the first place. Combining these three facts, we conclude 
that at present there is little evidence that modifications 
in the atomistic structure of the experimental contacts 
could exist, which are not properly included in the the- 
oretical modeling, and therefore would hint to an expla- 
nation as to why the theoretical conductance exceeds the 
experimentally measured one so much. 
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FIG. 7. Transmission of benzene- 1,4-dithiol calculated 
within DFT (solid Au54, as in Fig. 6; dot-dashed Aul4), 
HF (dashed Aul4) and extended Hiickel (dotted Aul4) with 
14 and 54 atom gold clusters simulating the electrodes. DFT- 
and HF-calcuations show large differences which indicates a 
breakdown of the EXCA explicated in section II B. 



Instead, we advocate another resolution of the puz- 
zle, namely that the mechanism described in section II 
becomes active because deviations from KS-states and 
HF-states are large. To illustrate that the proposed 
mechanism may indeed induce quantitative changes by 
orders of magnitude, we compare the transmission calcu- 
lated within the approximations DFT, HF and extended 
Hiickel. Fig. 7 exhibits the expected behavior: the 
DFT-transmission peaks are shifted as compared to the 
Hartree-Fock result and overlap much more strongly. In 
addition, the DFT-transmission has a broad background 
contribution only weakly depending on c;nergy. It stems 
from electronic states that reside in the leads and "leak" 
into the molecule. The DFT-calculation overestimates 
tunneling of these states as well. The net result is a differ- 
ence in the zero bias conductance by two orders of mag- 
nitude depending on whether the DFT-approximation or 
the HF-approximation is being used. Note, that because 
of the theoretical argument given in section II B it is not 
entirely clear, apriori, that it is necessarily DFT that 
gives the better approximation. 

In order to reduce the numerical effort we have re- 
stricted ourselves to small electrodes consisting of 14 Au- 
atoms each. By comparing the DFT-results for Aul4 and 
Au54 depicted in Fig. 6 one can convince oneself, that 
this simplification leaves unaffected our main conclusion. 



IV. CONCLUSIONS 

The numerical resiilts presented in the proceeding sec- 
tion are in full agreement with expectations based on the 
theoretical analysis performed in section II: while the 



conductance of a Au-chain agrees well with experimen- 
tal results, large quantitative discrepancies for organic 
molecule attached to gold electrodes via a thiol-bond ex- 
ist. Arguments have been presented to the extent that 
these discrepancies cannot plausibly be explained by in- 
sufficient modeling. 

The standard DFT-approach can be exact only in the 
case, where the ground state is represented by a sin- 
gle Slater determinant which in turn implies that HF 
is exact. For the molecule this is not the case, however. 
Since an approximation for the ground state based on 
the Slater-determinant of KS-orbitals is a mixture of the 
wavefunction of the exact ground state with its excita- 
tions, the character of the approximate state tends to 
be too delocalizcd. Consequently, the molecule-load cou- 
pling comes out too strong and the level broadening is 
overestimated. We have given analytical arguments in 
favor of this picture. 

In order to improve upon the standard DFT-approach 
to transport the exchange correlation potential should be 
replaced by the appropriate potential for the quasistatic 
nonequilibrium in order to include the corrections due to 
dynamical screening that are ignored otherwise. Unfor- 
tunately, at present such a potential is not known. How- 
ever, our work suggests an alternative approach which is 
based on the Kubo formula. It enables us to to calculate 
the dynamical polarization of the molecule together with 
parts of the leads. Thereof we can obtain the dc-current 
in the limit of zero frequency. Work in this direction is 
under way. 

Finally, let us emphasize that our results have implica- 
tions for the calculation of nonequilibrium effects^^ like 
the polarization of those electronic states that do not 
carry the current. Two effects that work in opposite di- 
rections occur. For a given current the voltage drop at 
the molecule is underestimated by DFT-based calcula- 
tions (since g is too large) and therefore so is the in- 
duced polarizing field. On the other hand, the charge 
response of localized states to this field is overestimated 
since they appear too metallic. Consequently, interpret- 
ing DFT based transport calculations for nonequilibrium 
effects is not straightforward and possible asymmetries 
in the I/V-characteristics could be blurred. ^^'^^ 
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V. APPENDIX 

We derive Eq. 4 from Eq. 5. The effective single parti- 
cle Hamiltonian of the scattering problem is of the form 
H = Ho H', where 
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(20) 



The Hamiltonian of the uncoupled left (right) lead has 
been denoted by /t^ {hn), the Hamiltonian of the cen- 
tral unit (molecule) is he and the matrices th,R are the 
couplings of the left and right leads to the central unit. 
If the molecule is uncoupled, the scattering states are of 
the form 



$i = (<^z,0,0) $r = (0,0,<^0. 



(21) 



When the coupling has been switched on, the new scat- 
tering states "ili^r can be obtained from a Lippmann- 
Schwinger equation 



(22) 



where E labels the energy of the state and the resolvent 
matrix Q{z) = [^; — H]~^ has been introduced. The lesser 
function is defined by (x=(x,t)) 



a<(x,x',i?) 



E 

a=l,r 



d{t-t') *„(x)*:(a;') /„ e 



■.E{t-t') 



(23) 



Using Eq. 22 a matrix representation can be derived 

g<(x,x',i?) = (x|[n-gH']g< [i + H'gt]|x') (24) 

with 



5l 
g< = I 

g< 



(25) 



and 



5^(x,x',i?) = ^ / d{t-t') U^Wiix') fi e'^(*-*'). 

(26) 

and a similar equation for g"^. In Eq. 3 only those ma- 
trix elements of S^are needed, that connect the left and 
right boundaries of the central unit. These elements are 
given by the term in Eq. 24 quadratic in the external 
coupling. Since 



H' go< H' 









(27) 
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with 



tng'^t^n and since Qo is diagonal 
one has for the matrix element representing the central 
unit 



(28) 



and thus recovers Eq. (4). A similar derivation has also 
been given by Brandbyge et al?"^ . 



1 J. Reichert et al, Phys. Rev. Lett. 88, 176804 (2002). 

^ M. A. Reed et al, Science 278, 252 (1997); C. Kergueris 

et al, Phys. Rev. B 59, 12505 (1999). 
3 X. D. Cui et al, Science 294, 871 (2001). 

J. Park et al. Nature 417, 722 (2002). 
^ W. Liang et al, Nature 417, 725 (2002). 
^ K. M. Roth et al, J. Vac. Sci. Technol. B 18, 2359 (2000). 
^ J. Chen, and M. A. Reed, Chem. Phys. 281, 127 (2002). 
* S. T. Pantclidcs, M. Di Ventra, and N. D. Lang, Ann. N.Y. 

Acad. Sci. 960, 177 (2002). 
^ E.G. Emberly, G. Kirczenow, Phys. Rev. Lett. 87, 269701 

(2001). 

1° J. Heurich et. al, Phys. Rev. Lett. 88, 256803 (2002). 
see also e. g. A. W. Ghosh and S. Datta, cond- 
mat/0303630vl 

Also, the KS-potential is not known exactly which intro- 
duces an additional problem, the "self interaction problem" 
(see e. g. R. M. Drcizlcr, and E. K. U. Gross, Density Func- 
tional Theory, Springer, Heidelberg (1990)). It leads to an 
estimate for the HOMO-LUMO-gap too small which im- 
plies a behavior too metallic. However, the effect depends 
on the extension of the wavcfunctions and wc expect that it 
is not of primary importance for molecules with extended 
TT-systems like those in Fig. 1. 
" Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

Y. Imry, Introduction to Mesoscopic Physics, Oxford Uni- 
versity Press (2002). 
1^ E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 



(1984); for a recent review see o. g. G. Onida, L. Reining, 
and A. Rubio, Rev. Mod. Phys. 74, 601 (2002). 

^® Gianluca Stefanucci, and Carl-Olof Almbladh, cond- 
mat/0310084v2. 

^'^ see e. g. G. D. Mahan, Many-Particle-Physics Kluwer Aca- 
demic Third edition (2000). 
D. S. Kosov, J. Chem. Phys, 119, 1 (2003). 

i^X. Gonze, and M. Scheffler, Phys. Rev. Lett. 82, 4416 
(1999). 

^° J. Reichert et al, cond-mat/0212272vl. 

P.S.Damle, A.W.Gosh, and S. Datta, Phys. Rev. B (Rap. 

Comm.) 64, 201403 (2001). 
22 M. Brandbyge et al, Phys. Rev. B 65, 165401 (2002). 
2^ Y. Xue, S. Datta, and M. A. Ratner, Chem. Phys. 281, 

151 (2002). 

O. Trcutlcr, R. Alihichs, J. Chem. Phys. 102, 346 (1995); 

K. Eichkorn et al, Chem. Phys. Lett. 240, 283 (1995). 
2S We use the BP86-functional; A. D. Becke, Phys. Rev. A 38, 

3098 (1988); J. P. Perdew, Phys. Rev. B 33, 8822 (1986). 
2® F. Evcrs, F. Woigcnd, and M. Koentopp, unpublished. 

K. Hansen et al, Appl. Phys. Lett. 77, 708 (2000); C. Un- 

tiedt et al, Phys. Rev. B 66, 085418 (2002); J. J. Palacios 

et al, cond-mat/0202375; Mads Brandbyge, Phys. Rev. B 

60, 17064 (1999). 
28 H. Sellers et. al, J. Am. Chem. Soc. 115, 9389 (1993). 

J. R. Reimers et. al, Ann. N.Y. Acad. Sci. 960, 100 (2002). 
^° H. B. Weber et al, Chem. Phys. 281, 113, (2002). 

Y. Xue and M. A. Ratner, cond-mat/0303179vl 
^2 K. Stokbro et al, Com. Mat. Sci. 27, 151 (2003). 

S.N. Yaliraki, M. Kemp, and M.A. Ratner, J. Am. Chem. 

Soc. 121, 3428 (1999). 

K. Stokbro, J. Taylor, and M. Brandbyge J. Am. Chem. 
Soc. 125, 3674 (2003). 
^® For a simple model system (quantum dot with on site inter- 
action plus leads) analytical calculations can be peformed 
and compared with DFT. Consistent with our assertions, 
large deviations in the strongly nonlinear regime have been 
found, see C. Romeike, Diploma Thesis RWTH Aachen 
(2002). 



9 



